Relaxation in quantum systems. Manifestation of the state-selective reactive decay 
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The efTect of state-selective reactive decay on the relaxation kinetics of quantum multistate sys- 
tems is studied in detail in the Bloch-Redfield approach (BRA). The results are applied to the 
analysis of this effect in radical pair recombination kinetics. The BRA is shown to be able to de- 
scribe quantitatively most important specific features of the recombination kinetics including those 
predicted by phenomenological treatment and by recently proposed approaches based on quantum 
measurement theories. 
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I. INTRODUCTION 

Recently some attention has been attracted to the dis- 
cussion of theoretical approaches for describing the ef- 
fect of state-selective reaction on the reactive quantum 
systems. The most popular example of such processes is 
spin-selective condensed phase physicochemical reactions 
of paramagnetic particles.^ The interest in this problem 
is inspired by new ideas, known in the quantum measure- 
ment theory,— i^ which can probably be useful in the anal- 
ysis of relaxation in quantum systems4^ At first glance, 
these ideas really provide some new insight into the spe- 
cific features of kinetics of state selective processes. In 
support of this conclusion recently some results have been 
presentedj^^^ which seem to demonstrate some new pos- 
sibilities of the new methods for treating experimental 
results. 

The discussion of these methods stems from some 
seeming limitations of the traditional approach looking 
rather heuristic. The main point of this approach is con- 
sidered to be the description of spin-selective reactivityii 
The description is based on the simple model for the re- 
activity supermatrix K operating on the spin density ma- 
trix of the system p. In this model, as applied to radical 
pair (RP) recombination with the rates ks and kt in the 
singlet {S) and triplet (T) states, respectively, 

kp^\{Ks[Ps.p]+ + ^T[PT,p\ + ] (1.1) 

where Ps and Pt are the operators of projection on the S 
and T states, respectively, and [P, p]+ = Pp + pP- The 
analysis based on the quantum measurement approach 
shows^-^ that the expression (|l.ip leads to the experi- 
mentally distinguishable underestimation the dephasing 
effect of reactivity. 

In this short paper we will analyze the state-selective 
reactivity effects on the relaxation in quantum sys- 
tems within the (traditional) Bloch-Redfield approach 
(BRA) 4"— The BRA is the basic approach in the mag- 
netic relaxation theory, but, in principle, it is very gen- 
eral and under quite realistic assumptions can properly 
describe a large variety of quantum relaxation processes. 



The results of our work demonstrate that formula of 
type of eq. (|l.ip is directly follows from the BRA. It is, 
first, derived for the simple three-state system (THSS), 
in which relaxation in the two state system is induced by 
reactive transitions to the third state, and then obtained 
[just in the form of eq. (|l.ip ] as applied to the case of 
spin selective RP recombination, by treating the reaction 
as spin-dependent transitions from excited (electronic) S 
and T states of the RP to the ground state. 

It is also shown that the BRA enables one to essentially 
generalize eq. (|l.ip by taking into account possible ad- 
ditional ST dephasing, accompanying the reaction. The 
generalized formula appears to be able to describe all 
predictions of the quantum measurement approach, dis- 
cussed in refs. [4,5], under some assumptions on relative 
value of reaction and dephasing rates. The effect of the 
interradical-distance dependent reactivity and exchange 
interaction is also briefly analyzed within the approach 
based on the stochastic Liouville equatio n^'^'^"^ which al- 
lows one to rigorously describe the manifestation of rel- 
ative diffusion of radicals. 



II. GENERAL FORMULAS 

Here we briefly discuss the main points of the BRA 
which are important for our further analysis of relaxation 
kinetics liSi^ For convenience, in our discussion we will 
use the frequency units for energies (i.e. h = 1) and the 
notation /3 = l/iksT). 

In the BRA the evolution of the system is described by 
the Liouville equation for the density matrix /)(i)>i^iii 

p = -i[Hs,p] + Rp, with [Hs,p\ = Hsp - pHs. (2.1) 

in which Hg is the non-fluctuating part of the Hamilto- 
nian H of the system and R is the relaxation supermatrix 
(the matrix in the Liouville space of matrix elements of 
the density matrix). Thus, the problem of description of 
the relaxation kinetics reduces to the evaluation of the 
matrix R and subsequent solution of eq. (|2.1[) . 



For the system in a bath the Hamiltonia n^°i^^ 

H = H, + H,i + Hi, (2.2) 

where Hi is the Hamihonian of the bath whose specific 
form is not important for our further analysis. The term 



Hsi=y $„(q)A„ 



(2.3) 



represents the interaction of the system with the bath, 
in which A„ are the Hermitian operators in the space of 
the system {\v), v = 0,1,2) and ^^(q) are the lattice 
coordinate dependent Hermitian operators of fluctuating 
amplitudes of interaction (see below). 

In what follows it will be convenient to represent op- 
erators A„ as sums of those AJ^, defined byiSil^ 



A„(t)=e*^=*A„e-'^=*=^ A 



r iuj^t 



(2.4) 



which describe transitions between different pairs of 
states of Hs [A„ = A„(t = 0) = Er^n]- 

In the sum (|2.4p the frequencies w^ are the differences 
of eigenstates Uj of the Hamiltonian Hg corresponding to 
eigenstates \j) (i.e. Hg\j) = iy.j\j)): 



and the operators AJ^ are defined by 

A; = (j|A„|/)<5^,;,„^^,. 



(2.5) 



(2.6) 



The operators and frequencies in eq. (12. 4p satisfy the 
relations (A^)'l' = A"'' and oj,^'' = —uj^ (see examples 
below) .iaa2 

In the BRA the relaxation rates are determined by the 
pair correlation functionsi^ii^ 

K+^,{t) = {'^n{m^'), if-„,(i) = ($„$„'(*)), (2.7) 

in which 

$„(i) = e'-f^'*$„e-*^'* (2.8) 

and the averages (...) are evaluated over the equilibrium 
lattice density matrix pi : 

{A} = Tt{piA) with pi = e-^"'/TT{e''^"'). (2.9) 

In addition we will introduce their Fourier transforms 



JL'i^) 



dtKt,{t)e^- 



(2.10) 



for which we get /+ ,(a.) = e^"J-,Ju;) = e^"J+,J-uj). 
The most important BRA results are conveniently ex- 
pressed in terms of Fourier transforms 

J„„'(^) = WnrA^) + Jn'.M\ (2-11) 

of symmetrized functions ivr„„/ (i) = 5[^r!n'(*)+-^r7'n(^)]- 



In particular, the supermatrix R is given byiSilS 

R=\{Ri+R2), (2.12) 

where 



R^p =^J„,„(a;;)r[A;,,p],A 



(2.13) 



R2P =^ J„'„(w;)tanh(i/3<,)[A«, [K',p\+\ (2-14) 

nn'r 

with [A, p\ = kp- pA and [A, p]^ = Ap + pA. 

Strictly speaking, in addition to relaxation the BRA 
also predicts some frequency shifts, ^"^^ though, we are 
not going to analyze them in this work. 

Concluding the short review of the BRA we recall the 
conditions of validity of this approach. The most impor- 
tant is the condition Tc <C ||i?|j— "— in which Tc is the 
characteristic correlation time of Hsi fluctuations, i.e. the 
average time of decay of correlation functions K^^,{t). 
This condition insures the validity of the second order 
perturbation approximation used in the derivation of the 
kinetic equation (|2.ip . Some additional condition is also 
required which ensures fast thermalization of the lattice, 
and in particular, fast relaxation of the non-diagonal di- 
agonal elements between s and I systems generated by the 
interaction Hsii^ This additional condition is, in reality, 
not very restrictivei^ and we are not going to discuss it 
in this work. 



III. THREE-STATE SYSTEM 

In this Section we will demonstrate the effect of reac- 
tivity on the quantum evolution kinetics by the exam- 
ple of the relaxation in the simple THSS, schematically 
shown in Fig. la. The THSS consists of the subsystem 
of two states (|0), |1)) and the third state |2), transition 
to which model the reaction. In the analysis we will also 
take into account the additional 0—1 dephasing [because 
of possible fluctuations of 0—1 splitting (see Sec. III.C.)], 
but, first, we will neglect this effect. 

In the considered THSS the terms of the total Hamil- 
tonian H [eq. (|2.2p ] are defined as follows 



H, = c^o|0)(0| + ic^,(|l)(l|-|2)(2|). 



Hsi = Hsi{12) 



2" 

^ ■'■ — 'n=x,y 



$n(q)o-„ 



(3.1) 

(3.2) 



In the operator H^i, representing (fiuctuating) interac- 
tion between |1) and |2) states, in which an, {n — x, y), 
are the Pauh matrices in the space (|1}, |2)): ax = 
|1){2| + |2)(1| and ay = i(|2)(l| - |1)(2|) . In principle, 
Hsi can be treated as the interaction of a quasispin, asso- 
ciated with matrices ax,y, with the fluctuating magnetic 
field, whose components are $n(q), {n — x,y). 

Formulas p.l|) and p.2p present the most general form 
of the two-state Hanfiltonian in the absence fiuctuating 
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FIG. 1: The schematic picture of the levels and transitions 
between these levels for: (a) the three-state system {wi and 
W2 are the rates of population relaxation, and woi is the — 
1 dephasing rate) and (b) radical pair {us and kt are the 
reaction rates, and kst is the ST^ dephasing rate, {fj, = 0, ±). 



diagonal interaction. As for the lattice coordinate de- 
pendent (fluctuating) components $„ (q) , their particular 
forms are not of real importance for our further analysis. 
The operators AJ^ introduced in formula (12. 4p are writ- 
ten as 

Al = zA,i =i(a,+za,) = i|l)(2|, (3.3) 



A^ 



-iA: 



i(a,-*a,)-i|2)(l| (3.4) 



As for fluctuations of $„, for simplicity, they are sug- 
gested to be uncorrelated and axially symmetric, i.e. 



Jn7i'{i^) = J{uj)Snn', {n,n' =x,y). 



(3.5) 



A. Relaxation kinetics in two-state (1, 2) subsystem 

First of all, it is clear that in the absence of interaction 
between the state |0) and two states |1), |2) the relaxation 
in the subsystem of two states |1), |2) is described by the 
conventional Bloch equations for the density matrix p in 
the subspace {|1), |2)}:^^ 

pil = -Wlipil + ■W22P22 (3.6) 

P22 = -W22P22 + Wiipii (3.7) 

Pm = -WnPM. (m = 12, 21). (3.8) 

The rates are given by 



Wii 



- p/3". 



eP^^W22^Jii^s)/il + e-P'^'), (3.9) 

Wn = Wi2 = ^{Wii+W22) = ^J{0Js)- (3-10) 

In particular, for large splitting ujq 3> l3~^ 

Wll =w^ = J{uJs), W22 = 0, Wn = ^Wll. (3.11) 



B. Relaxation kinetics in three-state system 

The analysis within the THSS shows that, in addition 
to relaxation in (|1), |2))-subspace, the fluctuating inter- 
action Hsi induces also the relaxation of — 1 and — 2 
phases, i.e. the decay of non-diagonal elements po^ ^nd 
Ppo with /z = 1,2. 

The expressions for corresponding phase relaxation can 
straightforwardly be derived with the use of eqs. (I2.12p 
and (Hm. 

For example, the decay rate of poi is determined by 
the matrix elements 

{0\Rip\l) - J(w,)(l|A°Jl)poi = -iJ(w,)poi, (3.12) 
and 

(0|i?2p|l) - J(c^,)tanh(i/3c.,)(l|AOjl)poi 

= -ij(ws)tanh(i;3ws)poi (3.13) 

where 

(3.14) 



A° =AiA-i + AiA-^ 

xy XX ' y y 



Similar formulas can be obtained for (l\Rjp\0) elements. 
Substitution of expressions (|3.12p and p.l3p into eq. 
(|2.12[) yields simple kinetic equations 

Pmm' = -WQiPfxfx', with p,p.' = 01, 10. (3.15) 

describing the exponential — 1 dephasing with the rate 

woi = ^w^i = y{uJs){l + e-f'^^yK (3.16) 

The analogous consideration enables one to derived the 
kinetic equations for relaxation of — 2 phases: 



PpA" == -wo2P^^', with /x^' = 02, 20, (3.17) 



where 



W02 = ^W22 ^ yiuJs){l + e^'^^yK 



(3.18) 



It is important to note that in the case of irreversible 
1 ^> 2 relaxation transitions, corresponding to the limit 
/3ujs -^ c», the obtained kinetic equations can be repre- 
sented in a simple matrix form with the use of the op- 
erator Pi = |1)(1| of projection onto the reactive state 
|1): 



p=->??[Pi,p] + , 



(3.19) 



with w^ deflned in eq. p. lip and [Pi, p]+ = Pip + pPi- 
It is easily seen that eq. p.l9p correctly describes relax- 
ation not only of the diagonal elements of the density ma- 
trix (population relaxation) but also of all non-diagonal 
elements (phase relaxation). 

Naturally, the effect of additional irreversible relax- 
ation transitions — ?► 2 with the rate w§g, can also 
be described, but with the more general equation p = 



-i(u;^[Pi,p]+ + u;o°g[Po,p]+), with Pq = |0)(0|. Of 
course this equation predicts additional dephasing in the 
system of type of that considered above. 

Note that similar expression, can be written for re- 
versible transitions (in the case of finite /3ws). One should 
only add the terms proportional to P2 = |2)(2| and the 
operator describing 1 — 2 (or — 2) relaxation. 



the relation Vi{p) — ^{PipQi + QipPi), where Qi — 
1-Pi. 

The operator (|3.24l) also predicts additional 1 — 2 de- 
phasing. This effect, however, is not relevant to the prob- 
lem under study and will not be discussed in this work. 



IV. RADICAL PAIR RECOMBINATION 



C. Effect of splitting fluctuations 

In the above analysis we have not taken into account 
possible fluctuating interactions diagonal in the bases 
(|0), |1), |2)) (the diagonal contributions to Hgi). 

In the two state system considered above (see Sec. 
IIIA) the effect of the fluctuating diagonal interaction 
of type of ^$z(q)(T^ is known to reduce to additional 
dephasing ,Sii£"i2 

In particular, in the simple case of isotropic fluctuat- 
ing interaction {^n}, {n — x,y,z), resulting in the in- 
teraction Hsl = ^Y.n=x,y,z^n{<i)(Jn, 1-6. for J„„' (w) = 

J{uj)5nn' , {n, n' = X, y, z), the dephasing rate w„ = W12 
is represented a o^^i^^ 



On = h{Wn + Wll), with 



J(0). 



(3.20) 



Comparison of this formula with eq. p.lO[) shows that 
the above-mentioned additional dephasing manifests it- 
self in the contribution ^w'^ — ^ J(0) to the dephasing 
rate. 

Naturally, similar effect is predicted in the presence of 
the additional fluctuating diagonal interaction, changing 
(0 - 1) splitting: 



i/.,(01) = i(|0)(0|-|l)(l|)a>oi(q) 



(3.21) 



Assuming uncorrelated fluctuations of $01 smd the field 
components $„ [see eq. p.2[) ] one obtains for the rate of 
0—1 dephasing 

WQi ^ ^(woi+wn), where woi = -'^Di(O). (3.22) 

The Fourier transformed correlation function Joi(w) can 
be calculated by substituting $01 (q) into eqs. (12. 7^ - 
(ETID . 

Similar expression can be derived for the rate of 0—2 
dephasing caused by the additional fluctuating diagonal 
interaction, changing (0 — 2) splitting. 

As in the absence splitting fluctuations, considered 
above, the obtained relaxation equations are conveniently 
combined into matrix equation which in the irreversible 
relaxation limit /3ujs — > cxd is expressed in the form 



where 



p=-UwrAPi,p]+ + woiAip)}. 



Vi{p)^\[Pi,p]+-PipPi. 



(3.23) 



(3.24) 



is the operator (in the Lindblat forni^) describing 0—1 and 
1 — 2 dephasing. This fact becomes especially clear from 



In this Section we will apply the results, obtained 
above, to the analysis of the irreversible-reaction effect 
on the kinetics of RP recombination. This process is 
known to be spin selective, i.e. the reaction rate depends 
on the total electron spin S = Sa + Sb of two radicals 
a and h. The corresponding states and transitions (as 
well as additional ST dephasing with the rate kst) are 
schematically shown in Fig. lb. 

The spin/space evolution of the RP is described by the 
spin density matrix p in the space of states of the total 
electron spin: the triplet states jPf), |T-), |Po) and the 
singlet state \S). In the simplest kinetic approach this 
density matrix satisfies the Liouville equatiortiii^ 



P= -i[Hp,p\ - Kp 



(4.1) 



in which Hp = Ha + Hi, is the spin Hamiltonian of 
the pair of non-interacting radicals (its particular form 
is not important for our further discussion) and K is 
the reaction/relaxation supermatrix, i.e the matrix in 
the space of matrix elements liyiy') = \v){i''\, {v,v' — 
S, Tfi with /i = ±, 0). 

In the BRA the general expression for K can be found 
by straightforward generalization of results, obtained in 
Sec. Ill [see eqs. p.23p and (|3.24j . In what follows we 
will consider the case of irreversible reaction in both \S) 
and |T^) states with rates ks and kt, respectively, and 
additional ST dephasing with the rate kst, expected to 
be induced by the fluctuating exchange interaction. The 
rate kx of reaction in [P^) states, as well as the rate kst 
of STfj^ dephasing, are assumed to be independent of the 
spin projection, i.e. of /u, (/Lt = 0, ±). 

Under these simplifying assumptions, in accordance 
with results of Sec. Ill, the supermatrix K can, in gen- 
eral, be written as 



Kp = ^{ks[Ps,p]+ 
= ^{ktp+ {ks 



kt[Pt,p]+ + kstVs{p)} (4.2) 
kt)[Ps,p]+ + kst'Ps{p)}, (4.3) 



where 



Ps = \S){S\, PT = l-Ps=Y.^^^jTf^)(T^\' (4-4) 
are operators of projection on S and T states, and 
Vsip) = ^[Ps,p]+-PspPs = \{PspPt + PtpPs) (4.5) 

is the term describing additional ST dephasing. Note 
that Ps(p) = Vrip) = \[Pt,p]+- PtpPt- 



In accordance with general analysis in Sec. Ill and 
formula (|4.2p the reaction/relaxation rates 



k„^, = {vv'\K\w'), [v,v' = 5,ro,±) 

are written as 

kss = Ks, kr^T^^KT, (Ai = 0,±) 
ksT,, = kT,,s = hii^s + kt + kst)- 



(4.6) 



(4.7) 
(4.8) 



Note that kr^r^^ and ksT^ are independent of /i. 

In our derivation of these expressions we did not spec- 
ify the final states for reactive transitions from S and T 
states. These final states are, quite probably, different 
for the initial S and T states. However, the choice of 
final states does not, evidently, affect the results in the 
considered case of irreversible reactions. 

The rates us, kt, and kst can, in principle, be eval- 
uated by means of the expressions, presented in Sees. II 
and III, in terms of correlation functions of correspond- 
ing interactions. Unfortunately, for many systems, of 
type of those discussed in our work, the obtained general 
expressions appear to be not very useful for accurate cal- 
culations of the rates because of the complexity of the 
systems which results, in particular, in the complexity of 
evaluating correlation functions K^^,. Nevertheless, the 
expressions are quite helpful for qualitative and semi- 
quantitative analysis, as it will be demonstrated below. 

It is easily seen that in the absence of additional ST 
dephasing {kst = 0) equation (|4.2|) reduces to the con- 
ventional formula (|1.1[) . thus showing that this formula 
is quite rigorous and is valid within the region of appli- 
cability of the BRA. 

As for the applicability of the BRA as applied to the 
RP recombination, according to the general comments 
presented in the end of Sec. II the BRA is valid when 
||iC||Tc <C 1, where Tc is the characteristic correlation 
time of fluctuations of Hsi . The time Tj, can roughly es- 
timated as the inverse characteristic frequency of vibra- 
tional and/or rotational motions which determine Hsi 
fluctuations. The most reasonable estimation for this 
time is Tc ^ 10^^'^s^^. Taking into account that usu- 
ally Ks, kt, Kst ^ lO^^s"^, we arrive at the estimation 
||iir||Tc < 10~^ which ensures good accuracy of the BRA. 
In fact, this inequality justifies quite natural require- 
ment that the probability of relaxation transitions to the 
ground state during the characteristic time Tc should be 
small enough. 

Of certain interest is the special property of the model 
(|l.ip . in which the density matrix p{t) [obeying eq. (|4.ip ] 
is represented in the form p{t) — \ip{t)){ip{t)\, with 
the "wave" functions \ip{t)) and {ip{t)\ satisfying the 
Shrodinger-like equations d\ip)/dt = —[iHp -|- \{nsPs + 
ktPt)M a.ndd{^\/dt = -(V'IH^p + iI^s-Ps + ktPt)]. 

It is worth noting that the limit of negligibly weak ad- 
ditional ST dephasing, when kst "^ ks,kt, can hardly 
be considered as quite realistic, because of expected small 
values of reaction rates in the case of large splitting uj^ of 



excited (reactive) S, T states and the ground state, which 
is typical for the considered process. The fact is that, 
according to eqs. p.23p and p.24|) . the reactive transi- 
tion rates k^, {v — S, T) , are determined by the Fourier 
transforms J^{ujs) of the correlation function of the inter- 
action, inducing corresponding transitions; Ki, ^ Jy{ujs)- 
The functions Jy{ujs) rapidly decrease when increasing 
iOsTc, where Tc is the correlation time of interaction fluc- 
tuations. Typically lOs is of order of the splitting of elec- 
tronic terms, while r^ (for vibronic coupling) is expected 
to be of order of (or larger than) inverse vibrational fre- 
quencies iOy <C Ws, SO that oJsTc ^ 1 and therefore we get 
for the total ST dephasing rate kst [see eq. (|4.8p ] the 
relation kst > ks, kt and therefore ksT,^ > kss, kT^T^- 
To illustrate these general arguments, note that typ- 



ically the rate ks 



W 



10 s and kt < ks- As 



for the dephasing rate kst, it can be estimated assum- 



ing by the simplest BRA relation*'"iA kst ^ J Tc, where 
J is the amplitude of fluctuations of the exchange in- 
teraction due to relative motion of radicals in the well 
of the attractive interradical interaction potential and 
Tc ~ 10~^^s. Assuming that J ~ Jo{aX), where Jq ~ 
10^^s~^ is the exchange interaction between radicals at 
a contact interradical distance r = d, a ~ 10®cm~^ is 
the rate of decrease of this interaction with the distance 
r: Je{r) = Joe-°'^'-'^\ and A > 10-^° cm is the esti- 
mated amplitude of stochastic relative motion of radicals 
at short distances. Substitution of all these parameters 
into formula for kst yields kst ^ 10"'^"'^ s~^ (> ks, kt). 

The additional arguments in favor of fast ST dephas- 
ing, corresponding to the relation ksT^, > kss,kT^T^, 
can be obtained by the analysis of the effect of stochas- 
tic (diffusive) relative motion of radicals with the use of 
the stochastic Liouville equationJ^ii^ The fact is that 
so far in our analysis we have not specified the mech- 
anism of stochastic motion along lattice q-coordinates 
which lead to fluctuations of H^i, though in our above 
estimations we have implied that the coordinates q cor- 
respond to the localized rotational/vibrational motion of 
the RP complex (cage) at the short distances r ^ d. In 
reality, however, very often one should take into account 
the additional coordinates describing relative diffusion of 
particles in the area of large space. In many case the 
contribution of these coordinates to the relaxation pro- 
cess in the RP appears to be more important than that 
of rotational/vibrational motion at short distances J^ii^ 

Naturally, diffusion leads to the stochastic change of 
the interradical distance r and therefore to fluctuations 
of distance dependent interactions and kinetic parame- 
ters. In describing the manifestation of relative diffu- 
sive motion the effect of reaction and relaxation at short 
distances can be modeled by reaction/relaxation term 
K{r)p similar to that in eq. (|4.ip but with the distance 
dependent supermatrix K{r), i.e. with the supermatrix 
(|4.2p in which the rates K^{r), {u = S", T, ST) are distance 
dependent. This model is expected to be fairly accurate 
since the characteristic scale of these r-dependences is 
a~^ ~ 10~®cm, while the typical jump length Id of the 



diffusive motion of not very small molecules is smaller 



than 



'-D ^ 10 cm. The diffusive relative mo- 



tion also leads to fluctuations of the exchange interaction 
J ex (r) at distances larger than the contact distance d (the 
distance of closest approach). The fluctuations of Jex{r) 
result, in turn, in those of ST splitting and therefore in 
ST dephasing. 

The analysis with the stochastic Liouville equation 
shows^ii^ that the efficiency of (distance dependent) 
reactivity and exchange interaction is characterized by 
the corresponding radii. In particular, the reaction radii 
Iviy (v — S,T), for S and T states resulting from reactiv- 
ities K^{r) — K°Ao(5(r — d), {1/ = S,T), localized in the 
narrow region ~ Ao near the distance of closest approach 
d, can be evaluated byi^ 



li^i^ = d[q^/{l + Qi,)], {v ^ S,T), 



(4.9) 



where qi, — n^{d\o/D) and D is the coefficient of rela- 
tive diffusion of radicals. As for the dephasing radius, it 
can be evaluated assuming exponential distance depen- 
dence Jeir) = Joe~°'^^~'^\ In the most realistic limit of 
relatively strong interaction, when Jq > Da^, the ST 
dephasing radius 1st is independent of dephasing at a 
contact and is given byi^ 

1st = lTs^d + a-i[1.14 + ln(2| Jol/Da^)]. (4.IO) 

i.e. for |Jo| > Da'^ the dephasing radius 1st > hs:lTT- 
Noteworthy is that since the rates of reaction in the 
states T^, {p = 0, ±) and dephasing rates, predicted by 
the supermatrix (j4.2|) . are independent of the projection 
H and therefore the reaction radii Wf^T,^ and Ist^ = It^^s 
are also independent of fj,: It^t^^ = Itt and IsT,^ = 1st- 

Most clearly the meaning of this relation between 
radii can be demonstrated as applied to diffusion as- 
sisted processes in cages: potential wells, micelles, etc. 
The kinetics of these processes is known to be close to 
exponential™"™ and the time evolution of spin density 
matrix is described by the kinetic equation of type of eq. 
(|4.ip with (diagonal) supermatrix K whose elements /c^^/ 
can be estimated in terms of corresponding radii^^'^^ 

{i^i^' \k\iyiy') = K,, « Dl,,,/Z, [v = S, To,±), (4.11) 



^drr^e '^"('') is the partition function of 



where Z = f 

particles within the cage [or in the potential well u{r)]. 
Recall that the radii It^t^, (m = 0,±), and IsT^^ are in 
dependent of /i so that fcr^T-^ and ksT^ do not depend on 
^ either: Ut^t^ = kTT, ksT^ = kT^s = ksT- 

Taking into account the above-mentioned relation be- 
tween the radii, we get ksT > kss,kTT- In terms of 
the rates of the THSS relaxation supermatrix (|4.2p . (14. 3p 
this inequality means that kst > ^^Sy^T- In reality, 
the dephasing rate ksT can essentially larger than kss 
and kTT since the corresponding radius 1st increases 
with the decrease of the diffusion coefficient D (though 
this dependence is very slow). For example, for realistic 
value&l^ Jo = 10^^ s~^ (low but still quite realistic es- 
timation), a = lO^cm"^, and D — 10^''^cm^/s one gets 



IsT — d « 4.14-10~®cm. Taking into account that usually 
dw 3-^5-10^^ cm we conclude that the rate of dephasing 
can be about twice as large as that of reaction. 

It is important to emphasize that in the realistic strong 
exchange interaction limit |Jo| > Da^ the radius 1st 
and, therefore, the rate ksT is independent of kst- This 
means that in the strong interaction limit the above dis- 
cussion of specific features of the superoperator K and, in 
particular, the relation between reaction and relaxation 
rates is not quite relevant. The fact is that in this limit 
the effect of dcphsing is determined by fluctuations of the 
exchange interaction at large distances r ~ 1st governed 
by relative diffusion of radicals. 

Note that different possible models of the recombina- 
tion process, including the models based on quantum 
measurement formalism, correspond to different relations 
between ksT , kss, and fcy t, and therefore can inter- 
preted within BRA type approaches under certain as- 
sumptions on values of interactions in the system. 

In the end of this Section it is also worth noting 
that the above discussion of the form of the reac- 
tion/relaxation supermatrix K is, of course, important 
for principle understanding the specific features of the 
relaxation kinetics of spin selective processes. As far as 
applications are concerned, however, the theoretical anal- 
ysis shows that for majority of processes in a wide region 
parameters the particular relation between elements of 
this supermatrix only weakly manifests itself in exper- 
imentally measured observables. For example, for the 
probability of the spin dependent and diffusion assisted 
RP recombination the effect of the specific relation be- 
tween Iss and 1st is characterized by the parameter 
£. = Q{lsT — IssY /D^— where Q is the characteris- 
tic spin dependent interaction in free radicals (hyperfine 
interaction, Ag change of the Zeeman frequences) The 
small value of this parameter means that the probability 
is insensitive to the particular relation between ^55 and 
1st- In particular, for realistic values D = 10~^cm^/s, 
1st — Iss — 3 • 10^^ cm, and Q — 10^s~^ we get fairly 
small ^ = 0.09, which means that at low viscosities (for 
highly mobile particles) the the effect of some difference 
between Iss and 1st can hardly be observed experimen- 
tally. As the viscosity increases, however, this effect be- 
comes quite measurable. 



V. DISCUSSION AND CONCLUSIONS 

The main goal of this short paper is to demonstrate 
that manifestations of state-selective reactivity in the ki- 
netics of relaxation of quantum systems can quite conve- 
niently and physically reasonably be treated within the 
BRA. It is important to note that the BRA is a math- 
ematically rigorous approach and absolutely transparent 
from physical point of view. The BRA allows one to 
analyze the contributions of different mechanisms and 
interactions to the kinetics of the processes giving deep 
insight into its specific features. 



The important example of such processes is spin- 
selective reactions of paramagnetic particles in which the 
reaction kinetics is essentially governed by (quantum) 
evolution of the spin subsystem via spin-dependent re- 
activity. 

In our works we have illustrated the obtained general 
results by the example of RP recombination. 

In the simplest assumptions on the mechanism of the 
spin-selective reactivity the BRA is shown to predict the 
conventional expression (jl.ll) for the reactivity super- 
matrix K . Its physical meaning can quite clearly and 
rigorously be demonstrated within the simple variant of 
the model (|3.ip . (|3.2p with the Hamiltonian Hs given by 
eq. dSH]) but Hsi == w(t)(|l)(2| + |2)(2|), where v[t) is a 
real fluctuating function with zeroth mean: {v(t)) = 0. 
Solution of the corresponding Schredinger equation for 
the wave function \ip{t)) — J2j=o i2'^jW\J) i'^ ^^'^ ^'^*-" 
ond order in v(t)t ^ 1 yields for the initial condition 
ai(0) = ao(0) = 1, 02(0) = 0: ai{t) « e-*"=*/2[l-Aa(i)], 
where Aa{t) = /J dti /q' di2w(ti)w(i2)e*"=(*i~*^\ and 
ao{t) — £"*""*. Thus one gets the following expres- 
sions for the density matrix elements averaged over v{t)- 
fluctuationsi^^ pii{t) = (|ai(i)p} « 1 - 2Re(Aa(i)) and 
poi(i) = (a5W«iW) « [1 - (Aa(t))]e*(-»-"=/2)*. Note 
that for relatively small t, satisfying the relation Tc < 
t ^ l/y^{v'^), where Tc is the correlation time of v{t)- 
fluctuations, we can write the relation 2ReAa(t) « wnt 
in which wn = 2/^ dt{v{t)v{0)) cos{uJst) is (1 -^ 2)- 
transition rate. The above formulas show also the rate of 
dephasing wqi, i.e. the rate of decay oi pQi{t), is given by 

Woi — ReAa(t)/t — hwn, in agreement with formula 

(ED. 

The analysis of the obtained simple expressions demon- 
strates that the relation wqi = i^n [i-^- eq. (jl.ip ] 
actually results from the evident fact that, if pii(t) ~ 
|oi(i)|2 - e-"""*, then poi(i) '- ai{t) - e-"'"*/2. 

The BRA analysis of general mechanisms, taking into 
account possible strong dephasing which can accompany 
the reaction process, has enabled us to essentially gener- 
alize formula (jl.ip . The generalized BRA permits natural 
physical interpretation of all results, obtained within re- 
cent "advanced" approaches ^^li^ in terms of clear phys- 
ical parameters and under reasonable assumptions. 

For example, recently proposed mechanism, based on 
some ideas of the quantum measurement theory, predicts 
the ST dephasing with the rate twice as large as that fol- 
lowing from eq. (|l.ip . i.e. with the rate kgj, = kss + krr 
rather than with fc^^^ — 2(^55 + kTr)-'^' This means, in 
particular, that in the case of only one reactive state, say 
singlet (S), ST phase is predicted to decay with the re- 
action rate: ksT = kss- Attempting to apply this result 
to real processes, one should keep in mind that, first, the 
applicability of the measurement theory^ to the processes 
at molecular scale is not quite evidenli^S and, second, any 
value of the ST dephasing rate can easily be interpreted 
in much more clear way by the analysis of fluctuating in- 
termolecular interactions within the BRA. Note also that 



just the relation ksr = kss is predicted by the most re- 
alistic diffusion theory, briefly discussed above. The fact 
is that according to eqs (|4.9|) - (|4.11|) in the limits of high 
reactivity ks{dX/D) 3> 1, and not very strong exchange 
interaction Jo/(-Da^) ^ 1 one gets kss = ks ~ ksT ~ d 
(in this estimation we also took into account that usually 
ad> 1). 

Similar comments can be added to the results of re- 
cent works^^ concerning the study of RP recombina- 
tion within the models, which can be considered as some 
variants of the spin boson model;^ usually applied in the 
theory of interaction of quantum dots. The model used 
in ref. [5] predicts that RP interaction results in the reac- 
tivity term of the form Kp ~ Vs{p) ~ [Ps, p]+ — '^PspPs 
which, as it follows from eq. (|4.5p . implies ST dephas- 
ing without reaction (recombination). This result means 
that in the model, used in ref. [5], the process reduces 
to ST splitting fluctuations, naturally, resulting only in 
ST dephasing. The author treats this dephasing as a 
set of measurement events which, in accordance with 
the Von Neumann "reduction postulate" (or "collapse 
postulate" ),^°'^S are associated with the instant loss of 
phase. This interpretation, however, does not look quite 
convincing. Another spin-boson type model (with an- 
other Hamiltonian) is considered in the work [19]. This 
model predicts for K formula (|l.ip which means (in terms 
of our general analysis) that in the corresponding fluctu- 
ating part of the Hamiltonian the diagonal part, describ- 
ing the fluctuations of splitting of terms, is neglected. 
In addition to the mentioned limitations of the analysis 
in works [5,6], note that as applied to the RP recombi- 
nation any spin boson type models, including those ap- 
plied in refs. [5,19] are of fairly restricted applicability 
(mentioned by the authors), for example, in the case of 
strongly unharmonic or stochastic (i.e diffusion like) rel- 
ative motion of particles. Note also that this kind of 
models is too complicated to be useful for the analysis 
of contributions of different kinds of interactions as well 
as the simultaneous effects of localized and delocalized 
motions. 

Moreover our analysis of the RP recombination within 
the BRA and diffusion approximation shows that this 
process can easily and rigorously be described by direct 
studying the kinetics of quantum relaxation without at- 
tracting any additional ideas whose applicability to small 
systems of molecular scale has not be rigorously demon- 
strated so far, and (what is most important) which, in 
fact, do not provide any new insight into the problem in 
addition to that obtained above from rigorously derived 
kinetic equations. 

Concluding this short discussion we would like empha- 
size once more that the effect of the state-selective re- 
active decay on relaxation kinetics of quantum system 
can fairly easily and accurately be estimated within the 
BRA. This estimation is usually sufficient for majority of 
interpretations. Of course, in some cases more detailed 
analysis within more specific models is required, but the 
models based on general ideas of the quantum measure- 



mcnt theory can hardly help in this kind of analysis. 
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